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Abstract of FR2827391 

The invention relates to a method of estimating a rainfall rate in the case of rain which only generates 
negligible attenuation using a bipolar radar. The inventive method comprises the following steps: the 
bipolar radar is Used to measure the differential phase PHI DP and the apparent reflectivity Ze along at 
least one of the polarisations H or V over a given segment [rO, rl] of a ray path in relation to said radai 
an estimate of value NO*, which is representative of the dimensional distribution of the rain drops, is 
determined using the difference in differential phase between rO and rl and using an integral of a 
fimction of the apparent reflectivity Ze along the length of segment [rO, rl]; the precipitation rate value 
is derived at a point using No* and the apparent reflectivity at said point. 
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TECHNIQUE POUR L'ESTIMATiON DE PLUIE AVEC RADAR LORSQUE UATTENUATION EST NEGLIGEABLE. 

(57) L'Invention concerne un proc6cl6 d'estimation d'un 
taux de precipitation dans le cas d'une plule ne generant 
qu'une attenuation negligeable, au moyen d'un radar bipo- 
laire, constltue par les etapes suivantes: 

- on mesure, au moyen dudit radar bipolaire, sur un In- 
tervalle [ro , r^] donn^ de rayon de trajet par rapport audit ra- 
dar, ta phase diff^rentlelie <t>np et la reflectivite apparente Zq 
selon au molns Tune des polarisations H ou V; 

- on determine une estimation de la valeur No' represen- 
tative de la distribution dimensionnelle des gouttes de pluie, 
k partir de la difference de la phase diff^rentlelle entre ro et 
r^ et a partir d'une Int6grale d'une fonction de la reflectivity 
ajpparente Zq le long de I'lntervalle [ro , r^]; 

- on deduit la valeur du taux de precipitation en un point 
a partir de Nq' et de la reflectivite apparente en ce point 
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La pr§sente invention est relative aux techniques meteorologiques 
5 pour lestimation d'un taux precipitant au moyen d*un radar. 

L'estimation du taux de pluie a partir de mesures de radar est 
generalement confrontee au probleme de I'attenuation de I'onde radar, et a 
celui de la variabilite naturelie de la pluie. 

Le concept d'un radar coherent a diversite de polarisation associe a 
10 un algorithme appele ZPHI a ete decrit dans WO 99/38 028, comme solution 
apportee a ces deux obstacles dans des conditions operationnelles. 

ZPHI est un algorithnne profileur qui, a la base, utilise en entree un 
profil de r6flectivit§ mesuree Za et une contrainte donn6e par une difference 
de phase differentielle Odp entre deux points r^ et d'une ligne de vis6e. 
15 A partir de ces mesures, on determine Tattenuation specifique A et un 

parametre connu sous le nom de No*, param^tre quantitatif de la distribution 
des tallies de gouttes. 

Le taux de pluie R ^ estimer est obtenu comme une fonction de ces 
deux param§tres, pour les deux raisons suivantes : 
20 - Tattfenuatlon sp6cifique A n'est pas soumise aux effets 

d'attenuation. Son utilisation dans Testimation de R permet ainsi de 
r6gler le probleme de Tattenuation. 
- le parametre No* suffit pour decrire la variabilite naturelie de la 
pluie. 

25 Ualgorithme ZPHI s'applique aussi bien en bandes X et C, bandes de 

frequence sensibles a I'attenuation. 

On vise ici maintenant un precede etendu a partir du precede ZPHI, 
lorsque I'attenuation n'est pas importante. 

Pour ce cas particulier, il est possible de prescrire un formalisme base 
30 sur des elements de base, mais plus adapts et beaucoup plus direct. 

L'algorithme ZPHI peut^galement etre applique a la bande S, pour 
laquelle Tattenuation n'est pas importante. voire negligeable. En effet. son 
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interdt subsists dans la mesure ou il reproduit « en direct » la variabllite 
naturelle de la pluie. Le recours au parametre A n'est alors pas forcement 
n§cessaire. 

L'invention a done pour but de proposer un proc6d6 simple, fiable et 
5 efficace d'evaluation du taux precipitant dans le cas d'att^nuation faible ou 
negiigeable. 

L'invention propose pour cela un precede d'estimation d'un taux de 
precipitation dans le cas d'une pluie ne generant qu'une attenuation 
negiigeable, au moyen d'un radar bipolaire. constitue par les etapes 
10 suivantes : 

- on mesure au moyen dudit radar bipolaire, sur un intervalle [ro. ri] 
donne de rayon de trajet par rapport audit radar, la phase differentielle Odp et 
la reflectivite apparente Ze selon au moins I'une des polarisations H ou V ; 

- on determine une estimation de la valeur No representative de la 
15 distribution dimensionnelle des gouttes de pluie, d partir de la difference de la 

phase diff§rentielle entre ro et ri et a partir d'une integrate d'une fonction de 
la reflectivity apparente Ze le long de I'intervalle [ro, ril ; 

- on deduit la valeur du taux de precipitation repr^sentant en tout point 
S partir de No* et de la reflectivite apparente en ce point. 

20 On propose egalement selon l'invention un dispositif pour I'estimation 

d'un taux precipitant comportant un radar bipolaire, ainsi que des moyens de 
traitement, led it radar comportant des moyens pour la mesure de la phase 
differentielle et de la reflectivite selon au moins une des polarisations H ou V, 
caracterise en ce que les moyens de traitement mettent en oeuvre les 

25 differentes etapes de traitement du precede selon le paragraphe pr6c6dent. 

D'autres caract^ristiques, buts et avantages de l'invention apparattront 
^ la lecture de la description detaiilee qui va suivre. 

Comme pour le precede ZPHI, le precede decrit maintenant requiert 
uri radar ^ diversity de polarisation coherent. Les donn^es d'entr^e sent le 

30 profil de r6flectivit6 Zh (ou Zv), ainsi que la mesure de la phase diff6rentielle 
entre la vole H et la voie V. 
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On vise ici S determiner le parametre No* (ainsi que le taux de plule) 
dans les cas oCi I'attenuation est negligeable, par exemple cas de la pluie 
mesur6e par la bande S, cas de pluies suffisamment faibles pour les bandes 
C at X, ou encore cas de la glace. 
5 Comme dans le proc6d§ ZPHI, le present precede integre deux types 

de donnees : la reflectivite Zh (ou Zv.ou plus g6neralement Ze) et la phase 
diff6rentielle Oop. 

Dans la mesure ou I'attenuation est negligeable, la reflectivite 
mesuree peut etre directement utilisee dans I'estimation du taux de pluie R. 
10 moyennant une estimation du parametre No* qui est precis6ment I'objet du 
present proc6d6. 

Le proc6d6 repose sur la relation « universelle » qui lie la reflectivite 
equivalente Ze (mm^m'^) et le taux de phase differentielle KDp(°km*^) : 



15 



20 



K 



No 



(1) 



Ou a et b sont des coefficients specifies par le modele de diffusion, qui 
dependent du type de precipitations (pluie ou type de cristal de glace) et qui 
sont legerement fonction de la temperature. La detemiination du type de 
precipitations peut etre operee par une methode de classification du type de 
celle decrite par Straka, Zrnie et Ryzhkov (2000). 

Plus generalement, le modeie de diffusion d6finlt la relation : 



(2) 



En integrant (1) ou (2) entre les deux homes n et ra du segment 
d'integation (ri<r2). on obtient : 



(3) 



25 
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Soit: 



^0 = 



4>> 



(4) 



II ne reste plus alors qu'^ transporter cette estimee de No' dans la 
relation « univeriselle » llant R a Ze : 



10 
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(5) 



pour obtenir I'estimee de R. 

Dans le cas general ou la relation Kop-Ze s'ecrit sous la forme (2). (3) 



devient 



^ ds 



(6) 



0 / 



On ri'a alors pas de solution analytique pour No*, mais une solution 
num6rique peut §tre trouv6e par m6thode num§rique iterative en utilisant 
comme valeur de depart pour No*(« first guess ») la solution (4). 

Dans le cas de la glace, qui n'att6nue pratiquement pas I'onde radar 
quelle que soit la frequence utilisee (bandes X, C et S), les limites 
d'application ne sont pas fix6es. 

Dans le cas de la pluie, I'application du proc6d§ est contrainte par une 
attenuation hegligeable entre les points ri et r2 (r2>ri). Pius specifiquement, 
les conditions d'application qui conduisent a une erreur de 3 dB sur No sont 
(pour T = ICC, et une loi de distribution de tallies de goutles de type gamma 
avec un param^tre de forme egal a 2) : 

Pour la bande X : [ODp(r2) - <^Dp(ri)] = ^ 4° 

Pour la bande C : [ODP(r2) - ODp(ri)] = < 10° 

Pour la bande S : [ODP(r2) - <l>Dp(ri)] = ^ 64" 
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Les domaines d'application principaux de I'invention sont : 

les memes domaines d'application que dans le cas de 
I'algoritlime ZPHI, ^ savoir I'estimation de la pluie sur les 
bassins versants pour la surveillance des crues et la gestion 
des ressources en eau. Cette application est valable pour tout 
type de pluie pour la bande S. ou aux pluies suffisamment 
faibles pour les bandes C et X. 
I'estimation des precipitations en phase glace. 
L'invention presente egalement de nombreuses autres applications en 
meteorologie. 

[1] Differential propagation phase shift and rainfall rate estimation. 
Sachidananda M. et Zrni6 D.S., Radio Science, 21-2, p.235-247 
(1986). 

[2] Polarimetric method for ice water content detennination, Ryshkov 

A.V., Zrni6 D.S. et BA Gordon, J. Appl. Meteor., p125-134 (1998). 
[3] The rain profiling algorithm applied to polarimetric weather radar, 

Testud J.E. Le Bouar, E. Obligis et M. Ali-Mehenni, J. Atmos. Oceanic 

Technol., 17, p332-356 (2000). 
[4] Bulk hydrometeor classification and quantification using polarimetric 

radar data: Synthesis of relations, Straka J.M., D.S. Zrni6, A.V. 

Ryzhkov. J. Appl. Meteor.. 39, pi 341-1372 (2000). 
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REVENDICATIONS 



1. Precede d'estimation d'un taux de precipitation dans le cas d'une 
5 pluie ne g6nerant qu'une attenuation negligeabie, au moyen d'un radar 

bipolaire. constitue par les etapes suivantes : 

- on mesure, au moyen dudit radar bipolaire, sur un intervalle [ro, ri] 
donne de rayon de trajet par rapport audit radar, la phase diff§rentielle a>DP et 
la reflectivity apparente Ze selon au moins I'une des polarisations H ou V ; 

10 - on d6tennine une estimation de la valeur No* representative de la 

distribution dimensionnelle des gouttes de pluie, ^ partir de la difference de la 
phase differentielle entre ro et T^ et a partir d'une integrate d'une fonction de 
la reflectivite apparente Ze le long de I'intervalle [ro, h] ; 

* 

- on deduit la valeur du taux de precipitation en un point k partir de No 
15 et de la reflectivite apparente en ce point. 

2. Proc^de selon la revendication 1, caracteris6 en ce qu'on deduit No 
directement a partir de la difference de phase differentielle entre ro et ri et a 
partir de {'integrate de la r6flectivit§ apparente Ze elev6e a un exposant 
choisi. 

20 3 .Precede selon la revendication 2, caract§ris6 en ce que I'exposant 

choisi est un exposant b, qui verifie : 



0 J 



OU Kdp est te taux de variation de la phase differentielle le long du rayon, a et 
b §tant specifies par le type de precipitation considere et par la temperature. 

4. Procede selon la revendication 3. caracterise en ce que le type de 
precipitation est un type correspondant parmi te type pluie et le type cristal de 
glace. 
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5. Precede selon la revendication 1, caracteris§ en ce que la fonction 
de la reflectivite apparente qui est integree est d§pendante de No , et en ce 
que le proc§d6 Inclut le fait de calculer S plusieurs iterations cette fonction, 
avec una valeur No* detemnin6e k iteration pr6c6dente. 

6. Precede selon la revendication 5, caract6ris6 en ce que I'on calcule 
No* a partir de la relation : 



ds 



ou Odp est la phase differentielle et ou F est une fonction verifiant la 
relation : 



n: 



oCj Kdp est le taux de variation de la phase differentielle le long du rayon. 

7. Proc6d6 selon la revendication 5 ou la revendication 6, caracteris§ 
en ce que Ton utilise, d une premiere iteration, la valeur de No* donn6e par : 



^0 = 



<^iDp(^2)-<I>Df(^l) 



8. Proc§d6 selon I'une quelconque des revendications pr6c6dentes, 
caract6ris§ en ce que Ton determine le taux de pluie R S partir de No* et de la 
r6flectivite apparente Ze par la relation : 



R = cM-%'' 
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9. bispositif pour restimation d'un taux de precipitation comportant un 
radar bipolaire. ainsi que des moyens de traitement, ledit radar comportant 
des moyens pour la mesure de la phase diff6rentielle et de la r6flectivite 
selon au moins une des polarisations H ou V, caract§ris§ en ce que les 
5 moyens de traitement mettent en oeuvre les differentes etapes de traitement 
du precede selon I'une quelconque des revendications pr§c6dentes. 
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